Enhancement of dark matter capture by neutron stars in binary systems 
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We study the capture of dark matter particles by neutron stars in close binary systems. By 
performing a direct numerical simulation, we find that there is a sizeable amplification of the rate of 
dark matter capture by each of the companions. In case of the binary pulsar PSR J1906+0746 with 
the orbital period of 4 hours the amplification factor is ~ 3.5. This amplification can be attributed 
to the energy loss by dark matter particles resulting from their gravitational scattering off moving 
companions. 

PACS numbers: 95.35.+d 95.30.Cq 



I. INTRODUCTION 

If dark matter (DM) particles interact with nucleons 
they must accumulate in stars [l], [2| . Their subsequent 
annihilation or collapse into a black hole inside the star 
may produce observable effects whose non-detection may 
be used to impose constraints on dark matter properties. 

In ordinary stars these effects are difficult to detect. 
Given the direct constraints on the DM-to-nucleon cross 
section 0] , the DM capture rate in ordinary stars is too 
low to eventually form a black hole. Likewise, the heat 
produced in annihilations of DM particles is by many or- 
ders of magnitude smaller than that resulting from the 
nuclear reactions. Only if the annihilation of DM par- 
ticles goes into neutrinos, the ones produced in the Sun 
can be in principle observed by the neutrino detectors 

The accumulation of DM is more efficient in compact 
objects like neutron stars and white dwarfs. Unlike or- 
dinary stars, these have no internal heat sources and the 
heat produced by DM annihilations can, in principle, 
be detected (8l-[l4|. Another possible observable conse- 
quence of the dark matter accumulation in compact stars 
is its eventual collapse into a black hole [1, [l3 - [l7| . Both 
the DM annihilation and collapse depend crucially on 
the amount/rate of the DM captured by the star, so an 
increase in the capture rate is a potentially important 
factor which may result in stronger constraints on the 
DM properties. 

In order to be captured, a DM particle has to lose 
a part of its energy and become gravitationally bound 
to the star. This energy loss may result from the DM 
scattering off the star nucleons. For this mechanism to 
work, a DM particle has to cross the star surface. An- 
other mechanism is operative in binary systems where 
the gravitational field is time-dependent. In this case DM 
particles may lose their energy by means of the gravita- 
tional interaction without crossing the star surface. In 
the context of the Solar system this phenomenon, known 
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as the gravitational slingshot, is routinely used to ac- 
celerate/decelerate spacecraft of interplanetary missions. 
This mechanism affects larger number of particles than 
the scattering off the nucleons. Deceleration of DM parti- 
cles by moving companions of binary systems may, there- 
fore, increase the number of gravitationally bound DM 
particles and amplify their capture rate. 

In the Solar system the effect of the planets (notably, 
Jupiter) on the dark matter ca ptu re by the Sun has been 
quantitatively studied in Refs. |18l42d | and was shown to 
be small. The smallness is due to the combination of two 
factors: (i) the inverse process (acceleration of DM parti- 
cles by the gravitational slingshot) prevents a substantial 
accumulation of DM and (ii) the Jupiter velocity is much 
smaller than that of DM particles, so only a tiny part of 
the DM phase space is affected. 

Both arguments are not directly applicable to close 
binary systems involving neutron stars. While the re- 
acceleration of DM particles still prevents an infinite DM 
accumulation, this effect may be less important in view of 
a higher capture rate by neutron stars as compared to the 
Sun. On the other hand, the velocities of companions are 
comparable to those of DM particles, so that the whole 
phase space of DM is affected. Thus, an unsuppressed 
effect of the companion motion on the DM capture rate 
in such systems may be expected. 

In this paper we calculated numerically the amplifica- 
tion of DM capture rate in close binary systems involving 
neutron stars. As a prototype we considered the binary 
pulsar PSR J1906+0746 composed of two neutron stars 
of approximately 1.3 Mq with the period of orbital mo- 
tion of 4 hours. In this system the velocity of the com- 
panions is 270 km/s, which is comparable to a typical 
velocity of DM particles in the Galaxy. 

We have found that in this binary system each of the 
companions captures DM at a rate 3.5 times higher 
than if it were an isolated star. The enhancement can be 
attributed to the combination of two factors: the velocity 
of the star relative to the DM distribution, and the scat- 
tering off the moving companions, i.e., the gravitational 
slingshot. One can see that the first factor, calculable 
analytically, is smaller than one and for the above binary 
system equals 0.57. Thus, the amplification due to the 
gravitational slingshot alone is about ~ 6. We also found 
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that the amphfication factor decreases when the orbital 
velocity of the companions becomes much smaller than 
the typical velocity of the DM particles, in agreement 
with earlier studies. 

The paper is organized as follows. In Sect. we de- 
scribe the numerical procedure used to calculate the dark 
matter capture rate in a binary system. In Sect. IIIII we 
present our results and conclusions. 



II. NUMERICAL PROCEDURE 

Since the velocities of DM particles and neutron stars 
are comparable, for most of the particle trajectories the 
interaction with both companions has to be taken into 
account. No analytic approximation seems possible in 
this case, and particle trajectories has to be calculated 
numerically. 

In principle, the numerical calculation of the capture 
rate is straightforward: one injects particles according to 
the Maxwellian distribution, traces them until they cither 
cross the surface of the neutron star and get captured or 
leave to infinity, and determines the capture rate from the 
fraction of captured particles. In practice, a substantial 
gain in calculational time is achieved by treating some 
parts of this process analytically. 

Consider first the motion of the binary system itself. 
For simplicity, we assume that it is composed of two neu- 
tron stars of equal mass M moving on a circular orbit. 
From the Newton equations one finds 
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where R is the radius of the orbit and w is the frequency 
of the orbital motion of the neutron stars. The motion 
of a DM particle in the vicinity of this binary system 
is most naturally described in units where lengths are 
measured in units of R and times are measured in unit 
of l/w. In these units the velocity of the neutron stars 
is w = ujR = 1. Assuming the parameters of the binary 
pulsar PSR J1906-I-0746, in the above units the velocity 
of light is c = 1127aji?, while a typical velocity of DM 
particles in the Galactic halo is udm ^ 270km/s ~ ujR. 
Note that the radius of the neutron stars is much smaller 
than i?. 



i?NS ~ 12 km ~ 2 X lO^^i?. 



(2) 



Let us now turn to particle trajectories. Defining two 
distance scales -Rmax ^ R and Rq <^ R, we separate each 
trajectory in three different regions: (i) the exterior of the 
sphere r = i?max which we will refer to as the injection 
sphere (ii) the small spherical regions of the radius Rq 
around the neutron stars, and (iii) the interior of the 
injection sphere excluding the vicinities of the neutron 
stars. 

In the region (i) the binary system may be approxi- 
mated by a single central mass 2M, so the particle equa- 
tions of motion can be solved analytically. Given the 



asymptotic velocity v^c and the impact parameter p, the 
parameters of the particle trajectory at r = i?max are 
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where v is the particle velocity and -0 is the angle be- 
tween the velocity and the internal normal to the injec- 
tion sphere. We generate the asymptotic parameters v^o 
and p according to the distribution 
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The extra power of v^o appears here because we are in- 
terested in flux rather than the density of particles. 

For each generated Voa and p we then calculate the 
parameters of the trajectory v and ip at the boundary of 
the injection sphere r — i?max- These parameters should 
be supplemented with the position of the entry point (two 
parameters) and the orientation of the projection of the 
velocity in the plane tangent to the sphere. We choose 
these additional parameters randomly. The resulting set 
of parameters fixes the initial conditions for propagation 
of the particle trajectory inside region (iii). 

Inside the injection sphere the spherical symmetry is 
no longer a good approximation. In this region we propa- 
gate the particle trajectories numerically until they either 
escape outside (beyond 200i?) with positive energy or get 
within the distance Rq from one of the neutron stars. 

In the region (ii), i.e., close to one of the neutron stars, 
the gravitational force of the remote companion can be 
neglected, and the problem once again becomes tractable 
analytically. Taking the initial conditions from numerical 
simulations, the trajectory can be continued inside r = 
Rq and the point of the closest approach to the neutron 
star can be determined. It is expressed in terms of the 
particle energy per mass E and the angular momentum 
per mass J in the reference frame where the star is at 
rest. 
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where the units R = lu = I were used. In this way the 
fraction of particles that get closer than a certain distance 
to one of the neutron stars can be computed. 

Since the probability for a particle to actually cross 
the neutron star surface in numerical simulations is very 
small in view of the small value of R^s, eq. ([2]), it is not 
realistic to measure the fraction of such particles directly 
with a reasonable accuracy. The problem can be avoided 
by noting that in the gravitational potential of a single 
star the cross section for crossing a sphere of radius Rg 
is equal to 



a{Rs) = 7rRi[l 
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where Uas is the asymptotic velocity of a particle. For the 
actual parameters of the neutron star system the second 
term dominates when R. One thus expects that 

for a sufficiently small Rs the cross section of crossing the 
sphere r — Rs scales linearly with R^. Since the depen- 
dence on the velocity factorizes in this regime, the linear 
behavior with R^ is maintained after averaging over par- 
ticle velocities. Measuring the fraction of particles cross- 
ing r = Rs a.t different values of Rs small enough for the 
linear scaling, one may then extrapolate to the actual 
neutron star radius. 

We note in passing that eq. ^ also explains why the 
capture rate by a single neutron star decreases with the 
increasing relative velocity of the star and dark matter: 
with the first term inside the brackets on the r.h.s. ne- 
glected, the cross section is proportional to l/v"^. De- 
pending on the direction of the particle velocity the effect 
of the neutron star motion can have both signs. Upon av- 
eraging over the Maxwellian distribution the net result 
is the decrease of the capture rate with the increasing 
velocity of the neutron star. 



III. RESULTS AND CONCLUSIONS 

In our simulations we did not fix the absolute normal- 
ization of the distribution of injected particles. Instead, 
for each value of the parameters we performed two iden- 
tical simulations: one with moving neutron stars and one 
with the neutron stars fixed to their positions at a given 
moment of time. In the second case all the effects due to 
the star motion disappear, and the capture rate equals, 
to a good approximation, twice the capture rate of a sin- 
gle isolated neutron star. The amplification factor thus 
equals the ratio of the number of particles captured in 
these two simulations. Clearly, this ratio does not depend 
on the absolute normalization of particle distributions. 

In each of the two simulations we measured the number 
of particles that got closer than Rg to one of the neutron 
stars, as a function of Rs- The results are presented in 
Fig. [H The parameters of the simulation were as follows: 
^max = lOR, Ro = 10~^R, and the number of thrown 
particles was 10*. One can see that the expected linear 
dependence with Rs is well reproduced. Fitting both de- 
pendencies with a linear function and taking the ratio of 
the two coefficients we find that the amplification factor 
equals 3.5 ± 0.1. 

The amplification factor depends on the parameters of 
the binary system, as it certainly has to go to one for 
slowly moving stars. The neutron star motion (assum- 
ing stars of equal mass) is determined by their mass and 
rotation period. When the trajectories of DM particles 
are expressed in the units of R and l/w, the change of 
the binary system parameters translates into the change 
of asymptotic velocity of the DM particles expressed in 
terms of Rlj, and the value of the capture radius that cor- 
responds to the actual radius of the neutron star. As we 
have seen, the latter parameter is irrelevant. Thus, the 



only important parameter is the ratio of the asymptotic 
velocity of DM to that of the neutron stars. Equal values 
of this parameter imply equal amplification factors. 
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FIG. 1: The dependence of the number of "captured" parti- 
cles (i.e., those crossing the sphere r = Ra) on the capture 
radius Rs for moving (upper line) and static (lower line) neu- 
tron stars. The points represent numerical simulation with 
1-a statistical errors. The lines show fits with the linear de- 
pendence. The actual radius of the neutron star is shown by 
the vertical dashed line. The number of injected particles was 
10*. 

We have calculated the amplification factors for several 
values of the binary system periods, assuming the masses 
of the companions equal M = 1.3Mq. The results are 
presented in Table. ID The maximum amplification corre- 
sponds to the case when the velocities of the companions 
and the asymptotic velocity of dark matter particles are 
comparable, while when the neutron stars move much 
slower than the DM particles, the amplification factor 
approaches 1. The amplification factor also decreases 
for short periods when the neutron star velocity becomes 
larger than that of DM particles, because the suppression 
due to the motion of the neutron stars relative to dark 
matter (cf. eq. (jS])) is no more compensated by the gain 
in the phase space of DM affected by the gravitational 
slingshot effect. 



period 


amplification 


4h 


3.5 


8h 


4.3 


16h 


2.8 


32h 


1.5 



TABLE I: The dependence of the amplification factor on the 
period of the binary system. 

To summarize, we have shown that the motion of the 
neutron stars in close binaries leads to the amplification 
of the DM capture rate by a factor of up to 3-4 depending 
on the period of the binary system. Thus, binary systems 
are favorable place to put constraints on dark matter 
properties that follow from its capture in compact stars. 
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